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MECHANISMS  AND  MODELING  OF  WIND- INDUCED 
LOW-FREQUENCY  AMBIENT  SEA  NOISE 


INTRODUCTION 

During  the  investigation  of  the  characteristics  of  ambient  sea 
noise,  it  has  been  observed  that  certain  portions  of  the  noise  spectrum 
level  are  related  to  local  wind  conditions.  Knudsen  curves  were 
established  in  the  later  1940* s  by  grouping  the  noise  spectrum  level 
according  to  sea  state.1  A  generalized  wind  noise  spectrum  was  com¬ 
piled  in  Wenz’s  survey  paper  and  was  supported  by  the  results  of  many 
measurements  in  the  frequency  range  above  200  Hz.2  Below  200  Hz  it 
was  believed  that  ambient  sea  noise  is  dominated  by  nonwind  dependent 
sources,  such  as  shipping,  biological,  and  industrial  noises.  However, 
recent  analysis  of  ambient  noise  data  below  10  Hz  indicates  that  the 
noise  level  in  this  part  of  the  spectrum  also  has  a  considerable 
correlation  with  windspeed.3  This  particular  relationship  has  not  been 
fully  explored. 

There  are  many  hypotheses  being  proposed  to  explain  the  causes  of 
ambient  sea  noise  attributed  to  surface  wind  conditions.  The  noise¬ 
generating  mechanisms  derived  from  theory  can  generally  be  summed  up  in 
three  categories.  One  treats  wind  turbulence  above  the  ocean  sur¬ 
face,  4,5,6  The  second  category  deals  with  surface-wave  interaction. 4,7-11 
and  the  third  treats  surface  wave  and  turbulence  interaction. 12  Because 
the  theoretical  formulations  of  such  problems  are  based  on  assumptions 
that  have  little  supporting  experimental  data  (particularly  in  an  ocean 
environment),  these  theoretical  results  are  not  adequate  to  interpret 
the  measurement  noise  data  collected  at  various  conditions  in  the  actual 
ocean  environment.  Definite  conclusions  cannot  be  reached  because  there 
is  no  direct  correspondence  between  measurement  and  theory.  Noise  data 
below  10  Hz  are  especially  lacking.  Furthermore,  attempts  to  use  one 
type  of  mechanism  to  explain  the  entire  noise  spectrum  appear  to  be  an 
oversimplification . 

This  report  reviews  some  of  the  basic  sound-generating  mechanisms 
attributed  to  ocean  surface  disturbances  caused  by  the  wind.  The 
report  derives  the  ambient  noise  level  attributed  to  such  causes  and 
compares  the  theoretical  results  with  measurements  obtained  in  known 
ocean  environments.  The  study  mainly  concerns  itself  with  the  wind- 
induced  noise  spectrum  below  10  Hz  because  a  reliable  ambient  noise 
data  set  with  supporting  environment  information  has  been  reported 
recently  and  is  available.13  When  considering  only  the  low  frequency 
range,  some  of  the  formulation  of  noise-generating  mechanisms  can  be 
simplified  with  physical  argument. 
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The  analysis  discussed  in  this  report  extends  from  previous  work, 
but  the  statistical  nature  of  ambient  sea  noise  is  emphasized  instead 
of  deriving  a  deterministic  model.4'11  It  also  takes  into  considera¬ 
tion  the  ocean  bottom,  water  depth,  and  sensor  depth  so  that  the  pre¬ 
dicted  theoretical  noise  spectrum  level  will  closely  correspond  to  the 
measured  data  obtained  in  the  actual  ocean  environment.  Physical 
arguments  will  be  used  to  simplify  some  mathematical  operations  and  to 
understand  the  connection  between  the  noise-generating  mechanism  and 
the  noise  spectrum  level.  This  also  provides  an  understanding  of  the 
conclusions  from  this  study  and  suggests  how  far  the  result  will  remain 
valid  if  conditions  change  from  those  assumed  in  the  derivation. 

A  wind- induced  noise  model  will  be  suggested  as  a  result  of  this 
study  to  highlight  the  important  parameters  of  such  a  noise  source  and 
to  provide  a  basis  for  the  computation  of  noise  fields  in  computer 
modeling.  It  also  can  serve  as  a  reference  for  data  comparison  and  as 
a  guide  for  future  planning  in  ambient  noise  measurements. 


FUNDAMENTALS  AND  FORMULATION 


The  mathematical  description  of  the  acoustic  wave  is  based  on  the 
physics  of  fluid  mechanics.11*  To  identify  the  important  parameters 
linking  the  ocean  surface  disturbance  caused  by  wind  with  the  ambient 
sea  noise  level  detected  by  a  sensor  in  the  water  column,  the  fundamen¬ 
tal  laws  of  fluid  mechanics  are  reviewed  here.  Formulation  of  the  wind- 
induced  ambient  sea  noise  is  based  on  this  derivation. 


According  to  the  conservation  laws  of  mass  and  momentum  of  fluid 
mechanics,  the  relationships  among  the  density  of  fluid  pf,  its  particle 
velocity  u1,  pressure  pf,  and  the  viscosity  v  (shear),  v*  (dilatation) 
can  be  described  by  an  equation  of  continuity. 


i£l  +  u»  i£l 

3t  i  ax. 


+  p' 


(i) 


and  a  force  equation. 


3u|  8u!  a 

pi  —  1  +  q  1  u ?  —  J*  =  _  _ 2 —  +  — 

at  j  3Xj  axj  axj 


-  +  [v1  d  .  +  2vd|  1 


where  i,j,k  are  coordinate  index  for  coordinates  1,2,3,  i.e., 
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Equations  (4a)  and  (5a)  describe  the  flow  mode  and  equations  (4b)  and 
(5b)  show  the  coupling  of  flow  mode  with  acoustic  mode.  Because  of 
equation  (4a),  flow  mode  is  assumed  irrotational  and  velocity  potential, 
0,  can  be  defined  as 


U. 

i 


-  ^ 


3x. 

i 


(6) 


The  equation  of  continuity  for  flow,  (4a),  becomes  Laplace  equation 


32  32  32  \ 

^7  *  3^2  *  TZ?)  *  =  ° 


(7) 


By  neglecting  the  effect  of  viscosity  in  the  flow  and  recalling  the 
irrotationality  assumption,  we  see  that  equation  (5a)  can  be  reduced  to 


po  at 


lu2 


♦  F(t)  , 


(8a) 


which  is  one  form  of  Bernoulli’s  equation.  For  acoustic  mode, 
equations  (4b)  and  (5b)  can  be  combined  into  a  form  as 


2  i  .2 

9.  P  .  i__ 

2  2  2  p< 
3x.  c  3t 

l 


32u.U. 


°U*  *  a  a  r  "l 

- L- L  _p  -  (u,u,)+  -  v’d  6.  ,  +  2vd.  . 

3x^3*j  o  3x^3x^  i  j  9xi9Xj  L  ij  ijj 


(8b) 

2 

where  the  compressibility  of  fluid,  p  /p,  is  assumed  to  be  a  constant, 
c2.  Later  in  this  report,  under  an  analysis  of  wind  turbulence,  equa¬ 
tion  (8b)  can  be  recognized  as  an  inhomogeneous  wave  equation  with  right- 
hand  side  terms  as  noise  sources  or  sinks.  For  example,  the  term 
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3x . 3x . 
i  J 


is  called  the  Reynolds  Stress  and  is  responsible  for  noise  generated  by 
turbulence.15  It  depends  on  the  order  of  magnitude  for  variables  in 
equation  (8b)  (initial  condition)  and  the  environment  of  the  disturbance 
(boundary  conditions).  Further  simplification  of  this  equation  can  be 
made  according  to  the  situation  where  it  applies.  The  above  derivation 
not  only  defines  the  notations  used  in  the  analysis  of  this  paper  but 
also  emphasizes  the  connection  between  flow  and  the  acoustic  wave. 


To  analyze  the  wind-induced  low  frequency  ambient  sea  noise  based 
on  the  fundamentals  of  fluid  mechanics,  a  simplified  scheme  of  ocean 
environment  (as  shown  in  figure  1)  will  be  used.  As  this  study  is 
mainly  concerned  with  the  low-frequency  acoustic  wave,  the  physical 
dimensions  will  be  scaled  with  reference  to  wavelength.  For  example, 
at  10  Hz  the  wavelength  of  the  acoustic  wave  will  be  150  m  whereas  the 
corresponding  wavelength  of  the  surface  wave  (the  gravity  wave,  the 
capillary  wave  will  not  be  appreciable  for  f  <  13.5  Hz)  is  about  0.017 
m  (A  =  2iTg/n2)  .  The  boundary  conditions  that  apply  to  flow  and  to  the 
acoustic  wave  can  be  approximated  in  a  different  way.  To  simplify  the 
mathematical  analysis,  the  whole  scheme  shown  in  figure  1  is  divided 
into  three  regions.  In  region  I,  the  direction  of  wind  is  assumed  to 
be  parallel  with  the  Xj,  x2  plane.  Since  6  is  larger  than  the  boundary 
layer  thickness  (6  >>  5.5  v/U  R  ),  the  wind  condition  above  6  can  be 
considered16  as  a  laminar  flow.  ^Region  II  contains  air-water  interface. 
The  depth  of  disturbance  in  this  region,  d,  is  much  smaller  than  the 
acoustic  wavelength,  but  is  about  the  same  order  as  the  wavelength  of 
the  surface  wave.  Region  III  is  the  water  column  with  a  semi - infini te 
bottom  layer.  In  figure  1,  p  and  c2  designate  the  density  and  compres¬ 
sibility,  respectively,  in  each  medium  with  subscript  "a"  for  air,  "b" 
for  bottom  layer,  and  no  subscript  for  water.  Channel  depth,  0,  and 
sensor  location,  s,  are  also  displayed  in  figure  1  so  the  computed 
noise  field  can  be  closely  compared  to  the  measurement  data. 


In  this  analysis,  the  behavior  of  fluid  in  region  I  can  be  simply 
described  by  equation  (4a),  and  the  wind  speed  condition  provides  a 
boundary  condition.  The  noise  is  mainly  generated  at  the  air-water 
interface  in  region  II.  The  induced  noise  occurring  at  region  II  will 
then  be  propagated  into  region  III.  Because  d  is  only  a  fraction  of 
acoustic  wavelength  in  water,  the  pressure  disturbance  on  the  air-water 
interface  can  be  regarded  as  the  boundary  condition  at  x3  =  d.  Since 
d  is  also  approximately  the  order  of  the  wavelength  for  the  surface 
wave,  the  hydrostatic  pressure  caused  by  the  undulation  of  the  ocean 
surface  will  be  negligible  at  this  depth.  The  derivation  that  relates 
disturbances  on  the  ocean  surface  to  ambient  sea  noise  is  based  on 
this  simplification. 


5 


Figure  1.  Analysis  of  Wind-Induced  Noise 
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NOI SE- GENE RAT I NG  MECHANISMS 
Wind  Turbulence 


The  existence  of  an  air-water  interface  will  cause  the  laminar 
flow  of  wind  to  change  to  a  turbulence  flow.  The  physical  processing 
of  this  phenomenon  is  rather  complicated  because  many  parameters,  such 
as  undulated  boundary,  temperature  difference  between  air  and  water, 
viscosity,  and  evaporation  rate  of  water  should  all  be  taken  into  con¬ 
sideration.17  However,  the  problem  we  are  concerned  with  here  is  that 
of  the  disturbance  occurring  in  the  1  to  10  Hz  frequency  range  where 
the  vorticity  size  is  small.  Some  understanding  about  isotropic  turbu¬ 
lence  can  be  utilized  to  simplify  the  mathematical  treatment  of  the 
wind  turbulence  problem  above  the  sea  surface.18 


Because  the  purpose  of  analyzing  wind  turbulence  is  to  determine 
the  pressure  distribution  at  the  interface,  equation  (8b)  can  be  used 
to  derive  the  relationship  of  pressure  and  fluctuation  of  velocity. 
The  acoustic  wavelength  is  very  large  compared  with  the  dimension  of 
the  turbulence  boundary  layer,  so  the  retardation  of  wave  will  be 
ignored  and  the  interface  will  be  considered  as  x3  =  0.  Therefore, 
equation  (8b)  can  be  rewritten  as 


3  u  .U  .  2 

-p  - -  p  — - -  (u  u  ) 

ao  3x.3x.  ao  3x.3X.  ai  aj 
i  J  x  J  J 


-P  s 
ao 


(9) 


By  assuming  3p/3x3  «  0  around  the  boundary  x3  =  0,  we  see  that  the 
pressure  distribution  can  be  expressed  as 


p(x1,x2>0,t) 


ao 

2tt 


/  /  / 

0  -  00 


scltidz 

|x-Z| 


(10) 


where  6  is  the  thickness  of  the  turbulent  flow.  Because  s  (Z,t)  is 
not  known  exactly,  only  the  statistical  property  of  p  can  be  derived. 
This  problem  has  been  treated  based  on  the  theory  of  isotropic 
turbulence. 19,20  The  following  are  the  major  steps  in  the  mathematical 
operation  and  are  adapted  with  present  notations  to  implement  physical 
interpretation  in  the  derivation  of  the  wind-induced  noise  level 
computation. 


First  assume  that  the  turbulence  field,  U  1 s,  is  isotropic 
(Kolmogoroff’s  hypothesis)  and,  then,  the  timeaverage  of  the  spatial- 
and-time  correlation  for  S  will  depend  only  on  the  relative 
separation18 
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and 


S(Z,t)  S'CZ'.f)  =  fCIz-z'-Uj?  |,  1  Z3_Z3l  »  7  lZ3+Z3^  »  (H) 

=  FtlZ-Z'-UjC.C.a]  , 

=  f[|y-Y'  +  |  .C * z]  , 

=  fC|y-y'  -  Ujd.e.z]  , 


where  yfs  is  the  horizontal  separation  between  the  boundary  point  and 
the  source  point  in  the  turbulence  boundary  layer,  U  is  the  mean 
velocity  of  the  flow  at  1/2  (Z3  +  Z$),  and  x  is  the  time  delay.  Then 
equation  (10)  can  be  extended  as  follows: 


R  (x.x'.t.t')  =  p~(x1,x2,0,t)  p[xj,x^,0,t') 


(12) 


p  2 

Hao 


--3-  // 


4tt 


S(Z,t)S(Z't‘) 
Ix-zl |x'-Z' I 


dZ  dZ' 


P  2 
ao 


f-  /  /  dZ3dZ5  /  / 
4  it  0  J  °  -  «> 


4tt 


dy  •  dY '  F[ Iy-Y'-UjCI »C,z] 


-,2  2  FT  -  72  '2 

JIx-yI  +  Z  JIx'-y'I  +Z- 


Pao2  ,6,  00  dy  dY '  fCy-Y'+x-x'-IKC.C.z] 

-T-  /  /  dZ  dZ'  /  /  - - 1 - 


J|y|2*z2J|y-|2  ♦  l'\ 


Applying  the  relationship  of  Hankel  transforms  yields 


/  d^dy'G  (y)G  (y')F(y-^-x)  =  4tt  /  KdKg  (K)g  (K)f(K)J  (K|x|),  (13) 
0  0  1  Z  0 


where  G  and  g  are  Hankel  transforms  pairs 
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g(K)  =  /  YdY  Jq(Ky)  G(y) 


G(y)  ■  /  KdK  J  (KY)  g(k)  . 
0  0 


Equation  (12)  has  the  form 


6  2  z  00 

R  CC.n.r)  =  2p  2  /  dz  /  d c  /  e'2Kzf (K.c, z) 

P  30  o  0  0  N 


Jo[K  Ju- uiT)2  +  n2]  , 


where,  based  on  the  assumption  that  the  disturbance  is  stationary  in  the 
Xj  and  x2  direction,  relations  of 

5  =  IVXJ  * 


n  =  |x2-x'|  , 


f(K)  =  /  YdY  J  (Ky)F(y)  » 


»  YdY  J  (Kr)  -Ka 
r _ o  _  e _ 

0  2  2  k 

Y  +a 


have  been  used. 


The  power  spectrum  of  p  can  be  derived  from  the  Fourier  transform 

of  R  : 

P 


frppU.n.o))  =  37  ]  RpCc.n.r) 


-1«T  J 

e  dt 


2p  *  «  ,  .  t  2z  °°  Av  -2KZ-,..  . 

*ao  I r  dz  -ib)C  /•  d  /  dKe  ^ffK.c.zj 

ir  i  U.  *  U.  {.  C  /-3 - ^ 

01  10  w/Ul  K  v  K2- (w/Uj) 2 
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•  cos  (n  K2-(w/U1)2  ) 
as 

/  J  [a(x2+b2)  ^2]cos  (xy)dx  =  (a2-y2)  ^2  cos[b  (a2-y2) 1//2]  for  o<y<a 
0  0 

Changing  the  variable 


and  replacing  f(K,£,z)  by  the  average  value  in  the  turbulence  boundary 

layer  f  ,  equation  (16)  can  be  further  simplified  to 
avg 


<j>  (£>n»w) 

PP 


4  P 


avg 


HU) 


6 

/ 


0 


-io)£/U 

zdze  1 


-M. 

e  U1 


1+X 


cos 


(17) 


where  6  is  the  displacement  of  the  turbulent  boundary  layer.  To  com¬ 
plete  the  integration,  the  mean  velocity  profile  of  is  needed.  White 
used  the  relationship 


U  6 


00 


(18) 


to  obtain  the  numerical  result  of  equation  (17).  It  agrees  well  with 
the  wind  tunnel  experiment .2 1 >22*2 3  For  the  problem  we  are  concerned 
with  here,  it  is  desirable  to  have  a  close-form  results  derived  from 
equation  (17)  so  that  the  noise  field  in  the  water  column  can  be  com¬ 
puted.  Since  it  is  not  easy  to  reduce  equation  (17)  to  such  a  form,  an 
alternative  approach  is  to  formulate  a  semi -empiri cal  expression  based 
on  its  numerical  result  and  the  measured  value  from  the  wind-tunnel 
experiment.20  The  final  result  is 


*pp  =  eao2  6u 


1.5*10 


r  *  ,  ~  2  U)6  3-j 

[  i+io  y-  ] 
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(19) 


-i (u)£/U  )  -  0.12u>£/U  -  0.55u>n/U  , 

c  c  «*> 

•  e 

where  U  is  the  convection  velocity.  Experimental  data  indicate20  that 

U  and  tree  flow  velocity  U  can  be  related  by 
c  ® 


U  -0.024  oj6/U 

-£  =  e 
U 


(20) 


Fourier  transform  of  equation  (19)  gives  the  wave  number  and  frequency 
spectrum: 


S  (a>,K) 
PP 


1.5*10'6p  2  u3<5 

...  _ ao  oo 


Ci+io'2 


00 


7T 

2 


0 . 12oj 
U 


0. 12a> 


U 


c 


Q  .55ai 
U 


+  K 


♦ 

1 


-i  r  0.55q?  2 
J  L  U 


(21) 


To  use  the  result  of  equation  (21)  in  an  ocean  environment,  know¬ 
ledge  of  wind  profiles  above  the  air-water  interface  is  needed  to 
determine  the  thickness  of  the  turbulent  boundary  layer.  In  general, 
for  an  aerodynamically  smooth  flow,  the  profile  is  found  by  experiment 
to  have  a  logarithm  form  above  the  turbulence  boundary  layer  (the  mean 
velocity  profile  in  the  turbulence  layer  was  assumed  to  follow  equation 
(18)  in  the  derivation  for  numerical  computation) . 6 #24  However,  for  a 
rough  surface  (such  as  the  ocean  surface  under  windy  conditions),  the 
thickness  of  the  turbulence  boundary  layer  probably  has  to  be  related 
to  wave  height  as  is  observed  in  a  recent  study.25  A  good  reference, 
therefore,  would  be  one- third  wave  height.  According  to  the  field  data 
collected  during  fully  developed  seas,  the  turbulence  layer  is  assumed 
to  be 

6-Hl/3  =  4-426  x  10‘3  U2’5  ,  (22) 


where  6  is  units  of  feet  and  U  is  units  of  knots. 26,27 

00 
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Whether  or  not  pressure  fluctuations  with  spectrum  of  equation  (21) 
cause  noise  generation  in  the  fluid  below  depends  on  the  wave  number 
vector  k.  Only  those  disturbances  on  the  surface  where  k2  <_  (w/c)2  will 
generate  an  unattenuated  sound  wave  into  the  water  column  of  region  III* 
Computation  of  the  noise  field  attributed  to  such  a  noise  generating 
mechanism  will  be  carried  out  in  a  later  section. 

Surface  Wave 


The  theory  of  the  mechanisms  for  wind-induced  ocean  surface  waves 
has  been  developed  extensively  by  Phillips  and  Miles.28*29  The  power 
spectrum  of  such  waves  modeled  by  Pierson  and  Moskowitz38  according  to 
the  observed  field  data  has  the  form 


I(«) 


2  -5  -B(u)  /u))^ 

ag  a)  e  o 


(23) 


where  - 

a  =  8.10-10  , 

6  «  0.74  , 

u  =  g/U  ,  and 
o 

U  =  wind  speed  reported  by  weather  ship. 

Some  of  the  basic  properties  of  the  surface  wave  can  be  explained 
by  the  linear  theory  as  derived  by  equations  (6)  and  (8a)  in  the  pre¬ 
vious  section.  If  only  the  gravity  force  is  involved  (gravity  wave), 
equation  (8a)  becomes 


*  n2  ~  n 

at  '  2  u  -  «x3  "  0  ’ 


(24) 


and,  taking  9/9t  from  it. 


L±  .  !_  cy 2  .  „  =  o  . 


at 


at  (v*r  +  8 


(25) 


For  small  disturbances,  the  effect  of  quadratic  term  9/9t  (V<f>)2  can  be 
neglected  and  2 

S  S  A  (b 

(26) 


14  +  &T~=  0  • 

17  9x3 
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By  means  of  equations  (7)  and  (26),  the  dispersive  relation  of  the 
gravity  wave  is  obtained,  i.e.. 


(27) 


which  has  a  velocity  potential  of 

,  .  lk,IVi(wt'VJ)  - 

<()(x1,x3,t)  »  /  b(u),k)  e  r  doidk  ,  (28) 

with 

|r|  =  Xj2+x22  and  |k12+k22|  =  k*  =  |kr2|; 


then  b  can  be  related  to  I(w)  by  matching  the  boundary  condition  at 
x3  =  C  (ocean  surface)  and 

|~  /  k_  b(oj,k)  dk |  2  =  I(w)  (29) 

Equation  (28)  also  indicates  that  pressure  due  to  flow  mode  (first- 
order  disturbance),  pQ  9<j>/8t,  will  be  attenuated  in  tbe  water  column 
for  x3  <  0.  Hence,  noise  detected  by  a  sensor  located  at  x3  <<  -l/k3 
is  not  caused  by  such  disturbances. 


The  second-order  effect  attributed  to  ocean  surface  disturbances 
can  be  derived  from  the  acoustic  mode  with  equation  (9)*  which  can  be 
rewritten  as 


dx. 

1 


c  3t/ 


-p 

o 


3x.3X. 

1  1 


(U.U  )  +  .  . . 

1  j 


(30) 


=  -p0  s(x,t)  , 


where  S(x,t)  is  the  source  term  caused  by  surface  disturbance  and  can 
be  derived  from  the  results  of  equation  (28).  Because  the  general 
solution  of  equation  (30)  has  the  form 


iwt-ik  *r  -  ik  *x 
_  .  _  r  3  3 

p  -►  e 


(31) 


•The  quadratic  term  in  flow  mode,  second  order  term,  is  put  back  in  the 
acoustic  mode. 
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with 


or 


2 

U) 

“T  = 


c 


» 


to  have  propagated  a  pressure  wave  in  the  x3  direction,  only  those  with 
|kr|2  <_  w2/c2  from  source  S(x,t)  excited  such  a  propagation  mode. 

The  wave  number  spectrum  of  the  ocean  surface  wave  is  not  totally 
understood  yet.  Experimental  data  and  proposed  theoretical  explana¬ 
tions  have  not  reached  satisfactory  conclusions.31’32*33  The  available 
empirical  model  is  given  by  an  angular  distribution,  such  as 


h(6)  =  2 


2s- 1  r 


l 3+1) 


2s 


r(2s+l) 


|  cos  -  (<f>-4io) 


(32) 


where  s  is  about  3  for  a  frequency  above  0.3  Hz  and  4>0  is 
of  the  wind.34  When  expressing  equation  (32)  in  terms  of 
Hfk)  becomes 


H(k) 


k  +  k-n. 


0.2 


9 


the  direction 
wave  number. 


(33) 


where  2 

|k I  =  w  /g 


and  is  a  unit  vector  along  the  axis,  or  the  direction  of  the  wind. 
Combining  equations  (33)  and  (23),  we  get  the  approximated  surface  wave 
spectrum  of  the  wave  number  and  the  frequency  is  modeled  as 


E(w,i()  =  i(w)  h(k)6(|k|  -  |k|)/|k| 


,  2  -5 

=  a  g  0) 


-0(a)  /(!)) 
0 


k  +  k 


_  3  (34) 

-^l)  6(|k|-|k'|)/|k|, 


where 

a'  =  1.62* 10‘3,  K'  =  w2/g  , 

and  other  notations  are  the  same  as  defined  previously. 
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Because  this  derivation  is  concerned  with  the  pressure  fluctuation 
in  region  II  with  a  length  scale  small  in  comparison  with  acoustic  wave¬ 
length,  the  equation  is  simplified  to 


S7  =  ^7  "W 


=  -P 


[£-  £-3  • 


0  3x.3X.  3x.  3x . 

i  J  i  J 


(35) 


Using  Fourier-Stieltjes  integral  representation  yields 


i(«t-k  • r)  +  ik.x 

p(x,t)  =  J  dP(k,u>)  e  r  J  5  dkrdio 


and 


4>(x,t)  =  /  d4(K,u)  e 


i(iut-K*r)  ♦  k.x 


3  3  ,r 


dkdu  , 


(36) 


where  r,  k,  K  are  considered  to  be  two-dimensional  variables  only;  i.e., 
r  =  (Xj»x  ),  etc.  Equation  (35)  then  reduces  to 


/  |k|  dP(oi,k)  e 


i(wt  -  k  *r) 


dkrdu 


■  -0  /  [k!2  K2  K.  ♦  K/  K.*  ♦  K.  K.  K.  K.  +  k!  K.  K.“] 

o  J  1  J  J  1  )  1  J  1  J  1  1  J  J 

(y'  \  ”, z"  "  iC(«,+»',)t-(K,*ic,,)r] 

•  d<P  (K  ioj  )  (K  )  e  v  v  7 


t»  r  2  M2  t  t  it  il  t  it  ft* 


(37) 


- »  -•*  *  it 

•  dK  dK  du>  du>  f 


by  evaluating  the  equation  at  2  0,  Let 


t 

w  ♦  Oi  *  <a  * 

„  t  .it 

K  +  K  *  kr  , 


(38) 
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Recall  that 


and 


|k|2  dP(«o,k)  =  -po  /  |V  (k.-K')  +  K.2  (kj“Kj) 2 

♦  h  Kj  °VV  <W  +  Ki  (ki‘Ki}  (kj'V2] 

•  d#  (K  ,w)  d4>'  (k-K*  .w-u/  )  dK*.  duj' 

■  ’o  1  [-kl  '‘l2  *  2k.  k2  Ki  *2  *  2i  kl  k3  S  S 

-  k2  *  2i  k2  k3  "2  *5  *  k3  Ki2]  d«V.») 


n  _  _  t  i  _  f  * 

d4>  (kr-K  , to— co  )  dK  du 


Ikl2  =4and  <2  +  ^2=  \j~\2  =  K 


'2 

3 


ik2ri<ik2i ,  ik2i«iK;2  +  k;2i  ikri2 « K;: 


(39) 


(40) 


and  consider  only  those  waves  that  propagate  downwards: 


I k3l > |kit  or  | k2 |  . 


Then  equation  (39)  can  be  further  reduced  to 

y}  f  3  t  f  t  „  f  ,  ,  f 

dP(a,k)~  p  — f-  /  | K  |  d$  (K  ,o)  )  d$  (k  -K  ,w-u>  )  dK  dw  .  (41) 
|k  |  *  r 

Because  of  the  condition  of  equation  (40),  it  may  be  assumed  that 


r  tt 

K  55  -K  , 

f  »» 

u>  ^  , 


(42) 


and 

with 


U)  =  2  00  =  2(u  , 


_  t 


dK  =  |  K  |  d|  K  |  d<f>  . 
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Equation  (41)  can  now  be  rewritten  as 


dPU.ic)  a  pocos20  j  K32  d$  (K  ,u>')d*  (-k'»u>')  6  (u)-2u/) dK *  du>' 


«)  dK* 


■  pocos20  /  K32  d»'(K  ,  y  u)  d$"(-K,y 

cos2e  /  tH* Ck’^1/2  [h"c-k')]  1/2  x  l4i'L;lib  dis' 


*  p 


-  P0  cos2e  (f)  i(|)  /  [h(*)3  '  [h (180°  ♦  ♦)]  '  d* 


1/2 


1/2 


=  po  cos2e  (|)  l(|)  0.8  /  sin 


tt/2 


.  3 


=  ¥p0  c-e  «(f) 


2  „/u>\ 


(43) 


where  the  relations 

2  2  2  lkJ 

(Kj)  (d$)  =*  0*0  I(w)H(K),  cos  9  =  *[Y| 

have  been  applied.  The  wavenumber  and  frequency  power  spectrum  at  the 
interface  is 

2 

r  04  2  42  4 

s  ((U,k)  =  ~~  p  cos  61  (u)ii>  ,  (44) 

pp  *5  o 

with  cos  9  = 

Equations  (43)  and  (44)  can  serve  for  the  boundary  conditions  ascribed 
to  the  surface-wave  noise-generating  mechanism  when  computing  the  noise 
field  for  region  III, 

Wave  Turbulence 


Wind  generates  waves  on  the  surface  of  the  ocean.  It  may  also  in¬ 
duce  flow  motion  at  the  surface  and  cause  vorticity.  Therefore,  there 
must  be  some  interaction  between  the  wave  motion  and  the  turbulent 
motion  at  the  surface  of  the  ocean.  Based  on  the  energy  balance 
between  amplification  of  vorticity  by  surface  flow  motion  and  dissipa¬ 
tion  of  vorticity  by  viscosity,  it  appears  that  wind-induced  turbulence 
is  small  (of  the  third  order).  8  However,  the  acoustic  wave  is  a 


17 


TR  5833 


second-order  quantity,  so  the  wind-induced  turbulence  may  make  some 
contribution  in  the  ambient  sea  noise  level,  particularly  if  it  inter 
acts  with  the  surface  wave  motion  that  is  a  first-order  quantity. 


To  evaluate  the  acoustic  pressure  due  to  the  interaction  of  wind 
induced  turbulence  and  surface  waves,  equation  (9)  is  rewritten  by 
including  the  relevant  terms: 


3 


2 


3xi3Xj 


*  -P 


o  3x.3x. 
i  1 


L3x. 

i 


(45) 


Applying  the  same  technique  used  for  surface  waves  by  the  Fourier- 
Stieltjes  integral  representation,  we  get 


p(x,t) 

4>(x,t) 

u. (x.t) 
J 

In  k,  w  space 

|  k  |  2  dP (<o,k) 


i (wt-k*  r)  -  ik  x  _ 
l  dP(k,w)  e  dk  du  , 

i (ut-K* r)  +  K.  x 

/  d4>(i(,u>)  e  4  dK  du  , 

i(tot-R*r)  -  iK  x 

/  dO  (K,o>)  e  *  dK  du  . 


-P 


o 


f  ? 

I  tv* 


»  t  I 

+K  K. 
j  i  3 


(k. 


i 


(46) 


1 2  1 
+  K.  (k.-KJ  + 
i  J  3 


_t  _  _  i  t 


Kicki'V  <yv]  "’(kV) 

,  (47) 

dK  do) 

dflCK*  JdU^  (k-ic*  jW-u*  )dK*  dj 


The  dispersive  relation  of  equation  (27)  does  not  apply  to  wind- 
induced  turbulence.  However,  simplification  of  equation  (47)  is 
accomplished  by  considering  that  there  is  no  correlation  between 
surface  waves  and  turbulence 
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k  \  4 

spp(u,,k)  a  po(tt)  /  Et(u-u)'  ,  k— K  * ) ,  d«\  die’ ,  (48) 


where 


»  .t  f  •  ® ( I K  |  -K  ) 

E  (u  ,K  )  -  I(u  )  H(K  )  - — 

w  Ik| 


f  .  .• 


Et(o)-(a)  ,k-K  )  =  dU^(o)-co  ,k-K  )dlK(u)-u>  ,k-K  ) 


- » 


the  spectrum  of  surface  turbulence,  according  to  the  measurement, 
generally  follows  Kolmogoroff 1 s  Law.35"37  From  the  field  data,37*38 
the  spectrum  of  turbulence  can  be  estimated  to  have  the  form 


E-,(a),K)  =  7  gy-=  ~ 

T  [l^2.104|Kr/3][l^] 


Using  the  results  from  equations  (23),  (33),  and  (49),  we  can  employ 
equation  (48)  to  compute  the  noise  power  spectrum  due  to  wave  turbu¬ 
lence  by  numerical  integration.  The  final  result  is 

-2  2  4  -2  lkJ 

S  (o>,k)  -  2.10  p  cos  e  w  u  cos  0  =  — —  .  (50) 

pp  iki 

NOISE  FIELD  COMPUTATION 

Disturbances  at  the  air-water  interface  cause  pressure  fluctua¬ 
tions  at  the  ocean-surface  boundary.  The  power  spectra  of  such 
pressure  fluctuation  attributed  to  three  different  types  of  noise¬ 
generating  mechanisms  have  been  derived  earlier.  The  noise  field  in 
the  water  column,  region  III  in  figure  1,  can  be  computed  according  to 
the  wave  equation, 

2  i  a2 

^Hr  ■  ~  f  =  0  »  (si) 

3x .  c  3t 


where  it  has  been  assumed  that  the  depth  of  surface  disturbance  d  is 
much  smaller  than  acoustic  wavelength  A.  Therefore  region  III  is 
source  free.  The  surface  disturbance  can  be  considered  to  be  a 
boundary  condition  imposed  at  x  8  0. 
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For  a  single  harmonic  disturbance  of  radian  frequency,  w, 
solution  for  equation  (51)  in  region  III  for  a  water  column  is 

ik  *r  -ik  x3  ik  x3 

o  t\r  ^  rm  »  3m  ,  A  3m  \ 
p  =  a  (k  J  e  (e  +  b  e  ) 

*m  m  m  m 


and,  for  a  bottom  layer,  is 


Pu  =  c  (k  )  e 
rbm  m  bm 


ikbrm*r  -  ikb3-X3 


where  ela)t  in  both  equations  is  the  suppressed,  r  =  (xj,X2),  k 
number  of  mode  m.  With  the  boundary  condition  at  =  0, 

ik  •?  ik  *r 

p  (k  )  e  srm  =  a  (k  )  e  rm  (i+b  ) 
sm  sm  mm  m 


and  it  results  in 


and 


at  X3  =  ° 


and 


where 


k  =  k 
srm  rm 


P  (k  )  =  a  (k  )(l+b  ) 
rsm  sm  mm  m 


ik  •  r 


a  (k  )  e 
m  m 


rm 


-ik  D  ik  D\ 

+  u  e  / 


=  cm  (kbJ  e 


ik  •  r  -  ik,  -  D 
b  rm  b  3m 


a  (k 
m  m 


ik  mr  (  -ik  D  ik  D^l 

)  e  rm  \-e  3in  +  b  e  3m  / 


i£,  -r  -  ik,  D 

=  -c  (k  )  e  brm  b3m  •  R  , 

m  bm 


R  = 


Pc 
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(52) 


(53) 

wave 

(54) 
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In  simplified  form 


k  =  kK 
m  bm 


(  “ik*  D 
le  3m  + 

(  -ik,  D 

\-e  3m  + 


ik,  D 

b  e  3m 
m 


■ik  D 
b3m 


*  -c  ft 
m 


-ik,,  D 
b3m 


and  equations  (54)  and  (55)  can  be  used  to  determine  a^  b  ,  and  c  . 
Hence,  the  pressure  field  in  the  water  column  at  a  particular  deptn,  s 
for  a  particular  wavenumber,  k^,  can  be  expressed  as 


p  (1+R)  e 


-ik  (S-D) 


+  (1-R)  e 


ik,  (S-D) 
3m 


cos  k3mD  ♦  R  sin  k^  D 


A  (k  )  OR2)  +  n-R2)  cos  2  k,  (S 


cos  k,  D  +  R  sin  k,  D 
3m  3m 

where  the  power  spectrum,  \s>  due  to  the  noise-generating  mechanisms 
caused  by  surface  disturbance,  has  already  been  derived  from  the  last 
section. 

In  general,  R  <  1  and  R2  <<  1 .  Equation  (57)  can  be  further 
simplified  to 

A  1  ♦  cos  2  k  (S-D) 

rm~  - 2 - ’  (58) 

cos  k,  D  ♦  R  S  k,  D 
3m  m  3m 
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The  total  field  intensity  at  frequency  u>  attributed  to  all  the  wave¬ 
number  contribution  will  be 


=  /  I  d  k 
m  m 


(59) 


Here,  since  the  noise  fields  of  different  wave  numbers  arriving  at  sen 
sor  location  s  originate  from  different  locations  at  the  surface,  inco 
herence  is  assumed  for  the  integration  expression  from  equation  (59) . 

The  dominant  contributing  wave  number  modes  in  equation  (49)  are 
those  that  meet  the  conditions 


or 


cos  k„  D  = 
3m 


k3m  =  (2n+1)  I  i 


(60) 


Those  modes  also  have  to  be  the  unattenuated  modes  generated  at  the 
surface.  They  must  satisfy 


2  -  (-)2  -  k2 

rm  \A  /  3m 


=  (“)2 

r,  x2n^)  ii 

\\  1 

L1  4  D  J 

or 

2D  1 

n  -  “  -  2  • 

(61) 


Hence,  the  intensity  of  the  noise  field  at  the  sensor’s  depth,  s,  in 
the  water  column  can  be  approximated  by  the  relation 


i  - 1 

n 


A(V 


1  ♦  cos  [(2n*l)  J  (1-  |)] 
_ 

R 


(62) 


where  n  must  meet  the  conditions  given  by  equations  (60)  and  (61) . 

The  corresponding  noise  power  spectrum  measured  by  the  sensors  at 
depth,  D,  will  be  r  *  c  1 

„  A(k.)  1  +  cos  (2n+l)  \  (1-  i) 

IW  -  ♦(-)  I  - 1  2  " - , 


(63) 
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where  A(k)  4>(w)is  the  noise  power  spectrum  generated  at  the  air-water 
interface  and  is  derived  in  the  last  section. 

MEASUREMENT  DATA  COMPARISON 

Figure  2  was  prepared  from  the  processed  low  frequency  ambient 
noise  data  collected  in  the  Bermuda  area  during  January  1966  and  illus¬ 
trates  the  data  set  reflecting  the  least  interference  from  a  nonwind- 
dependent  noise  source.13  The  plotted  noise  spectrum  level  was 
obtained  from  a  bottom-mounted  hydrophone,  900  m  (2600  ft)  below  the 
surface.  The  noise  data  are  grouped  by  corresponding  wind  speed  and 
show  distinct  profiles  above  1  Hz.  For  the  spectrum  level  below  1  Hz, 
the  data  becomes  questionable  because  of  system  noise  interference  and 
deterioration  of  system  response. 

The  theoretical  values  of  ambient  noise  attributed  to  three  hypo¬ 
thesized  mechanisms  are  computed  according  to  equation  (43)  with 
appropriate  parameters  (bottom  loss  data  are  about  8  dB) ,  and  the 
results  are  plotted  in  figure  3  for  comparison  with  the  measured  data 
at  20  to  40  knot  wind  speeds.39  Although  none  of  the  theoretical  values 
match  exactly  with  the  measured  data,  they  appear  to  be  of  the  same 
order  of  magnitude.  The  spectrum  shape  of  wind-induced  flow  turbulence 
noise  particularly  seems  to  follow  that  of  the  measured  data.  This  may 
indicate  that  some  of  the  constants  assumed  in  the  derivation  of  the 
generating  mechanism  are  not  quite  correct,  but  the  general  relation¬ 
ship  in  the  analysis  still  seems  to  hold.  It  appears,  judging  from  the 
shape  of  each  individual  curve  with  respect  to  frequency,  that  the 
wind  turbulence  and  the  surface  wave  are  the  dominant  cause  for  wind- 
induced  ambient  noise  below  7  Hz.  Wave  turbulence  becomes  an  important 
noise-generating  mechanism  for  frequencies  above  7  Hz. 

To  further  identify  the  important  parameters  of  the  noise-generat¬ 
ing  mechanisms,  figure  4  is  a  replot  of  all  the  measured  data  and 
theoretical  values  with  wind  speed  as  the  coordinate.  The  noise  level 
change  caused  by  the  surface  wave  is  small,  as  most  energy  in  wind  is 
converted  into  low-frequency  surface  waves  when  the  wind  speed 
increases  (Pi erson-Moskowitz  Model30).  On  the  other  hand,  the 
noise  level  change  caused  by  airflow  turbulence  and  turbulence  wave 
interaction  covers  a  wider  range  and  follows  almost  the  same  pattern  as 
the  measured  data.  At  5  Hz ,  these  curves  indicate  that  wind  turbulence 
and  surface  wave  noise  generation  mechanisms  provide  a  reasonably 
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corrected  value  of  ambient  noise  level.  At  10  Hz,  the  wave  turbulence 
noise-generation  mechanism  provides  very  good  predictions  on  the  noise 
level  in  comparisons  with  the  measured  data. 

Comparisons  between  theory  and  measurement  are  also  made  for  data 
collected  at  Grand  Bank,  as  shown  in  figure  5. 3  Agreement  of  the  order 
of  magnitude  still  appears  to  hold,  but  the  theoretical  values,  in 
general,  are  lower  than  the  measurement  results.  In  this  particular 
area,  it  seems  that  the  surface  wave  mechanism  is  the  dominant  noise 
source. 

MODELING 

Although  the  theoretical  value  based  on  the  analysis  is  not  com¬ 
pletely  in  agreement  with  the  measured  data,  the  order  of  magnitude 
and  spectrum  shape  provide  reasonable  support  for  the  hypothesis  that 
the  three  different  noise-generating  mechanisms  are  responsible  for 
wind-induced  ambient  noise.  The  main  discrepancy  is  caused  by  lack  of 
fundamental  information  on  air  turbulence,  wave  directionality,  and 
surface  wave  turbulence.  Experimental  data  for  these  parameters  in  the 
ocean  environment  are  sparse.  However,  a  semi-empirical  model  for 
noise  can  be  constructed  by  matching  the  spectrum  level  with  measured 
data  in  terms  of  the  parameters  identified  through  theoretical  analysis. 

There  are  other  variables  in  formulating  a  noise  model.  First, 
according  to  equation  (42),  bottom  loss  is  a  very  important  factor  in 
the  final  value  of  the  noise  field.  For  a  more  accurate  computation, 
a  computer  model  of  propagation  should  be  incorporated  in  the  noise- 
field  calculation.  Other  factors  affecting  the  noise  level  are  the 
geographic  location  of  the  site  and  the  time  history  of  wind  conditions. 
For  example,  if  the  measurement  site  is  close  to  the  coast  or  a  storm 
center,  it  is  highly  probable  that  waves  propagated  in  the  opposite 
direction  can  exist.  In  this  case,  the  dominant  mechanism  may  be 
caused  by  the  surface  waves.  Another  deviation  from  the  analysis  would 
be  a  short  duration  of  a  certain  wind  speed,  or  a  relatively  fast  wind 
speed  change.  The  fetch  will  then  be  small  and  the  stationary  nature 
of  the  surface  disturbance  cannot  be  assumed.  Probably,  the  predicted 
noise  spectrum  level  will  be  overestimated.  The  shift  in  wind  direc¬ 
tion  will  cause  a  change  in  the  wave-directional  spectrum  and,  hence, 
the  noise  spectrum  level. 

Formulation  of  a  noise  model  for  computer  operation  can  be  based 
only  on  the  mean  of  a  stochastic  process  of,  in  a  sense,  a  determinant 
model.  Since  the  means  are  well  established  for  computing  the  sound 
field  from  a  given  surface  source  distribution,  the  desired  noise  model 
is  the  noise-source  distribution  at  the  surface  (region  II  in  figure  1) . 
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From  the  above  consideration  and  the  results  of  the  noise-mechanism 
analysis,  the  suggested  form  of  surface-noise-distribution  model  can  be 
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where 


F(6)  =  noise  radiation  pattern  (cos  0),  and 

f,  g>  h,  ot,  6,  X  =s  constants  to  be  determined  empirically  that 
are  functions  of  wind  duration,  direction, 
and  site  topography. 

The  first  term  is  attributed  to  the  surface-wave-generating  mechanism. 
The  second  and  the  third  are  caused  by  wind  turbulence  and  wave  turbu¬ 
lence,  respectively.  More  data  are  needed  for  systematic  analysis  in 
order  to  validate  these  controlling  constants  of  the  noise  model. 


CONCLUSIONS 

Three  noise-generating  mechanisms  have  been  examined  at  the  air- 
water  interface  of  the  ocean  surface.  With  the  propagation  condition 
and  boundary  effect  incorporated,  the  derived  noise  field  based  on 
those  three  mechanisms  provide  a  reasonable  explanation  of  the  measured 
noise  data  for  various  wind  conditions.  The  physical  argument  used  in 
the  derivation  of  wind-induced  low-frequency  noise  enables  simplifi¬ 
cation  of  the  mathematical  operation  and  provides  a  better  under¬ 
standing  of  the  observed  phenomenon.  A  noise  model  can,  therefore,  be 
formulated  from  the  derived  relationship  with  appropriate  adjustment 
of  parameters. 

There  are  some  basic  problems  required  for  further  study.  The 
exact  turbulent  structure  above  and  below  the  surface  of  the  ocean  is 
not  fully  explored  for  various  wind  conditions.  Such  knowledge  is 
needed  to  extend  the  current  analysis  to  nonstationary  conditions  of 
ocean  surface;  i.e.,  the  transition  stage  of  a  passing  or  distant  storm. 

The  results  from  the  current  analyses  indicate  that  the  noise 
spectrum  level  is  sensitive  to  bottom  properties.  This  means  wind- 
induced  noise  might  be  geographically  dependent.  There  are  no  available 
data  measured  at  other  locations  to  verify  such  conclusions.  Some 

measured  low  frequency  data  have  been  published  but  there  is  a  lack  of 
supporting  wind  speed  time  series  information.40*  4i  The  Grand  Banks 
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data  show  a  different  slope  and  higher  noise  level  when  compared  with 
the  Bermuda  data.3  It  appears  that  the  surface-wave-generating  mech¬ 
anism  is  the  dominant  source  for  the  Grand  Banks  data.  Since  the  bottom 
loss  has  never  been  measured  at  either  site  for  the  1  to  10  Hz  fre¬ 
quency  range,  no  definite  conclusion  can  be  reached  on  its  effect  on 
the  measured  ambient-noise  level. 

It  has  also  been  hypothesized  that  wind-induced  deep  sea  currents 
might  be  another  noise-generating  source.  However,  the  observation  in 
the  Bermuda  area  indicates  that  the  magnitude  is  very  small  (lqss  than 
1  knot)  with  no  noticeable  wind  speed  correlation. 42  The  hydrophones 
used  for  noise  measurements  in  the  Bermuda  area  were  laid  or  placed  on 
the  bottom.  Therefore,  there  should  not  be  any  strumming  problem  with 
the  cable,  and  any  turbulence  around  the  hydrophone  would  probably  be 
small  with  such  a  small  current. 

Some  other  noise-generating  mechanisms  (such  as  microseismics, 
seaquakes,  shipping,  and  biological  effects)  may  also  be  responsible  for 
low-frequency  noise;  however,  they  are  not  correlated  with  wind  condi¬ 
tion.  The  analysis  discussed  here  does  not  include  such  a  study. 

The  result  of  the  current  study  may  be  extended  to  other  frequency 
ranges  of  the  noise  spectrum,  but  this  would  require  further  examina¬ 
tion.  Because  both  airflow  turbulence  and  surface-wave  generated  noise 
have  a  high  order  of  negative  slope,  their  contribution  to  noise  at  a 
high  frequency  may  be  small  and  surface  wave  turbulence  may  be  assumed 
to  be  the  dominant  mechanism.  But,  in  general,  from  10  to  200  Hz, 
other  non-wind-dependent  noise  sources  mask  the  whole  spectrum.  Wind- 
dependent  noise  becomes  appreciable  only  above  200  Hz.  It  has  been 
postulated  that  bubbles  and  water  drops  are  the  causes.43*44*45  Further 
study  certainly  is  needed  to  explore  such  a  possiblity. 

The  theoretical  analysis  discussed  in  this  report  does  not  result 
in  a  precise  picture  of  what  occurs  in  a  real  ocean;  however,  it  has 
identified  some  of  the  more  important  parameters  in  the  wind-induced 
noise-generating  mechanism  and  it  also  serves  as  a  guideline  for 
future  noise  measurement.  If  such  information  as  wind  direction, 
surface  turbulence,  wave  direction,  and  wave  turbulence  can  be  col¬ 
lected  simultaneously  when  making  noise  measurements  the  ocean  will 
be  better  understood. 
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